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This  document  is  a final  report  of  work  performed  under  Air  Force 
Office  of  Scientific  Research  Contract  No.  F49620-77-C-0121 . The 
contract  was  let  in  response  to  the  unsolicited  proposal  "Study  of 
Radiative  Scattering  Concept  of  Plume  Infrared  Radiation  Obscuration" 
submitted  to  the  Air  Force  by  Utah  State  University  in  March  1977.  The 
work  was  performed  between  1 August  1977  and  30  October  1978.  A three 
month  no-cost  extension  enlarged  the  original  1 year  contract  to  15 
months . 
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ABSTRACT 

Theoretical  and  experimental  studies  were  performed  to  measure  the 
synergistic  absorption  of  hot-gas  radiation  by  an  aerosol/absorbing  gas 
mixture  and  to  verify  quantitatively  the  infrared  absorption  of  carbon 
particles.  Computer  calculations  for  a pure  scatterer/NjO  absorbing  gas 
mixture  revealed  a synergistic  effect  of  only  10Z.  In  carbon  - NjO 
mixtures,  synergism  was  negligible.  Attenuation  measurements  were  made 
for  pure  carbon  and  TiO^  aerosols,  pure  NjO  samples,  and  the  combination 
of  each  aerosol  in  NjO.  In  all  cases,  the  net  transmittance  of  the 
gas/aerosol  was  described  by  the  simple  product  of  the  N^O  alone  and 
aerosol  alone  transmittances . No  synergism  was  detected.  Transmission 
calculations  for  carbon  show  it  a very  effective  candidate  for  plume 
obscuration . 
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I.  INTRODUCTION 

A combined  experimental  and  theoretical  (modeling)  program  to 
determine  the  effectiveness  of  aerosol  scattering  to  reduce  infrared 
radiation  from  jet  engine  exhausts  was  carried  out  under  AFOSR  contract 
No.  F49620-77-C-0121 . The  experimental  work  was  done  at  the  Aerospace 
Corporation  in  El  Segundo,  California  and  the  theoretical  effort  at  Utah 
State  University  in  Logan,  Utah. 

The  primary  infrared  emission  features  of  jet  engine  exhaust  plumes 
are  the  COj  and  H^O  vibration/rotation  bands  at  2.7pm  and  the  CO^ 
fundamental  band  at  4.3pm.  The  intensities  of  these  emissions  are  very 
high,  particularly  for  aircraft  operating  in  the  augmented-power  mode 
where  they  may  be  an  order  of  magnitude  higher  than  at  military  power 
levels,  and  provide  the  primary  signal  for  the  guidance  of  heat  seeking 
missies. 

The  natural  COj  and  content  of  the  atmosphere  can  provide  a 
significant  reduction  in  the  emission  intensities  over  long  atmospheric 
paths.  However,  because  the  emission  bands  arise  from  a 
high-temperature  environment  (T  = 900°K  for  military  power  and  up  to 
2100°K  in  the  augmented  power  mode)  they  are  substantially  wider  than 
the  corresponding  absorption  bands  of  the  relatively  cool  ambient 
atmosphere.  This  widening  is  caused  by  the  population  of  both  higher 
vibrational  levels  and  higher  rotational  levels  of  the  molecules  in  the 
plume  than  in  the  atmosphere.  Whereas  the  atmospheric  CC^  and  HjO  are 
very  effective  at  attenuating  the  plume  emission  band  centers, 
essentially  no  attenuation  is  provided  in  the  band  wings  and  this 
"wing-leakage"  radiation  can  be  propagated  over  many  miles  of  atmosphere 
with  little  reduction  in  intensity.  Over  very  long  atmospheric  paths, 
the  radiated  spectrum  approaches  a limiting  form  of  very  sharp  intensity 
"spikes"  at  the  band  edges,  particularly  on  the  short-wavelength  side  of 
the  bands.  The  principle  spectral  regions  of  concern  are  thus  the  red- 
and  blue-spike  regions  of  the  2.7  and  4.3pm  bands. 


I 


l 


I 


t 


( 


A previous  study  [Harris  et  al.,  1975]  considered  the  obscuration 
of  this  spike  radiation  by  gaseous  absorption.  In  concept,  a cool 
gaseous  shroud  would  be  formed  aerodynamically  around  the  plume.  The 
shroud  would  contain  a gas  or  mixture  of  gases  that  had  strong 
absorption  in  the  spike  spectral  regions.  An  experimental  search  was 
made  for  suitable  shroud  gases.  Several  candidates  were  found  which  had 
reasonable  obscuration  potential  for  one  or  more  of  the  four  spike 
features.  The  prime  examples  were  HjO  and  methyl  isothiocyanate  for  the 
NjO  4.6pm  CO 2 band  wing  and  NjO  and  formic  acid  for  the  2.9Pm  HjO/CC^ 
band  wing.  Estimates  of  the  amount  of  these  materials  ( hundreds  of 
pounds)  that  would  be  required  to  provide  a 90Z  reduction  of  intensity 
for  a period  of  about  10  seconds  indicated  that  the  concept  was 
marginally  feasible. 

Infrared  radiation  reduction  was  reported  for  aircraft  plumes  when 
an  aerosol  suspension  of  particles  was  formed  as  a shroud  around  the 
aircraft  exhaust  plume  or  in  the  plume  itself.  Obscuration  was  reported 
for  carbon  suspensions  [Varney,  1972;  Blozy,  1974;  Birstein,  1976],  oil 
mist  shrouds  [Varney,  1972],  and  purportedly  for  jet  fuel  mists  [Varney, 
1972;  Caldwell,  1975].  The  fact  that  these  techniques  employed  an 
aerosol  suspension  of  particulates  or  liquid  drops  would  suggest  that 
the  effects  of  radiation  scattering  play  a role  in  the  overall 
obscuration  effects  observed.  It  should  be  clearly  stated,  however, 
that  pure  scattering  phenomena  could  not  account  for  an  intensity 
reduction  when  the  field  of  view  of  the  observing  sensor  is  large  enough 
to  include  the  entire  plume  and  aerosol  shroud.  For  a cylindrical 
shrouding  geometry  in  vAiich  the  plume  is  surrounded  on  all  sides  by  the 
aerosol  shroud,  the  extinction  effects  of  side  and  backscatter  by  a 
particular  small  volume  of  the  shroud  must  be  cancelled  when  the 
scattering  contribution  of  all  such  scattering  volumes  is  considered. 

(A  small  percentage  of  radiation  could  conceivably  be  "funnelled"  to  the 
rear  of  the  cylindrical  .shroud  if  this  end  is  not  filled  in  by  turbulent 
closing  of  the  shroud.) 


The  only  real  mechanism  for  radiation  reduction  is  absorption.  In 
the  case  of  carbon  suspensions,  the  absorption  ia  evident  since  carbon 


is  black  throughout  the  infrared  region  of  interest.  For  the  oil  mist 
and  fuel  mist  techniques,  the  suppression  mechanism  is  not  entirely 
clear  since  neither  jet  fuel  (kerosene)  nor  presumably  the  oils  used 
provide  significant  absorption  in  the  infrared  regions  of  interest.  The 
suspected  mechanism  for  these  suppressants  is  that  they  do  provide  some 
absorption  and  that  the  effects  of  scattering  tend  to  enhance  this 
absorptive  potential. 

The  specific  question  addressed  in  the  present  work  is  whether  or 
not  the  superposition  of  a gas  and  aerosol  suspension  might  result  in  a 
synergistic  enhancement  of  absorption  such  that  the  net  transmittance 
through  the  mixture  is  less  than  the  product  of  the  transmittances  of 
the  gas  and  aerosol  taken  independently.  Both  theoretical  and 
experimental  efforts  were  undertaken  to  evaluate  any  synergistic 
absorptions  that  might  be  produced  by  scattering  particles.  Two 
particle  types  were  considered;  titanuim  dioxide  and  carbon.  The 
transparency  of  titanium  dioxide  in  the  infrared  allows  a study  of 
purely  scattering  interactions  in  an  absorbing  gas.  ^0  was  chosen  for 
this  gas  because  of  its  relatively  good  abscuration  in  both  the  2.9  and 
4.6um  wing  regions  and  because  of  its  ease  of  handling.  The  other 
particle  type,  carbon,  was  chosen  because  of  its  absorptive  as  well  as 
scattering  properties.  A brief  description  of  the  computer  code  used  to 
theoretically  calculate  scattering  is  followed  by  results  aimed  at 
evaluating  the  synergistic  effects  of  scattering  particles  and  absorbing 
gasses,  and  calculations  of  the  radiation  transmission  through  carbon. 
Results  of  radiation  measurements  through  scattering/absorbing  mediums 
are  presented  next  followed  by  the  conclusions  of  this  work. 


II.  THEORETICAL  EFFORT 


A.  Computer  Code 

The  theoretical  work  performed  under  this  contract  dealt  with  the 
calculation  of  infrared  radiation  through  a scattering  and/or  absorbing 
medium.  A computer  code  based  on  the  Braslau  - Dave  solar  radiation 
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model  [Braslau  and  Dave,  1973]  and  modified  by  Luther  [1974]  was 
obtained  from  the  University  of  California,  Lawrence  Livermore 
Laboratory.  This  code  was  developed  to  calculate  the  vertical  radiation 
transport  in  the  scattering/absorbing  portion  of  the  lower  atmosphere  of 
the  earth.  Modifications  were  made  at  Utah  State  University  to 
calculate  the  infrared  transmission  through  a one-dimensional  horizontal 
homogeneous  path  filled  with  scattering  and  absorbing  particles. 

The  basic  radiative  transfer  equation  for  monochromatic  radiation 
has  the  form 


dl( x ! p > $) 

j m ~ w(  ^ )J(  1 5 p,  $) 


(1) 


where  KxJp.^)  is  the  intensity  of  the  radiation  at  position  t in  the 
direction  p - cos  where  e is  the  angle  with  respect  to  the  normal  to 

the  slab  and  4 is  the  azimuth  angle  referred  to  an  arbitrary  meridian 
plane  (see  Figure  1).  Forward  (backward)  direction  is  represented  by  +4 
(~p).  The  optical  depth  t is  the  sum  of  optical  depths  due  to 
absorption  and  scattering  processes  and  is  defined  to  be  0 at  the 
source,  increasing  into  the  medium.  Values  for  the  scattering  depth  can 
be  computed  by  Mie  theory  [Dave,  1974].  Absorption  optical  depths 
depend  upon  the  gas  absorption  cross  section.  The  albedo  of  single 
scattering  w(t)  is  defined  as  w(t)  * Ax8/A-  the  fractional  change  in  t 
due  to  scattering. 

The  source  function  J(x;p,$)  is  given  by 
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J(  T » Jl»  T * M *♦} 


(2) 


where  P(T;y,4;ti'  ) is  the  phase  function  which  describes  the  amount  of 

radiation  scattered  into,  the  direction  (u,+)  from  the  direction  (p',41). 


Fig.  1.  Schematic  diagram  for  radiation  transport  calculation 


I 


,♦’)  - T(t)  M (v,$,u'  ,♦')  + (1  - T(  t)  ) R C M , ♦ ; u a ) (3) 

where  the  so  called  turbidity  factor  T(t)  is 
T(t)  - AT(s'm)/(AT(8'm)  + AT(s'r)) 

the  fraction  of  total  scattering  due  to  Mie  scattering.  The  quantities 
M and  R are  the  Mie  and  Rayleigh  scattering  phase  functions  respectively 
to  be  described  later.  The  B function  for  our  radiation  model  stands 
for  either  a unidirectional  source  or  an  isotropic  radiator.  A 
unidirectional  source  is  similar  to  solar  radiation  in  the  earth's 
atmosphere  and  can  be  written 


B(t;Mo**o)  “i  e T/M°  • (4) 

2 

where  irF  is  the  incident  flux  (watts/cm  -pm)  normal  to  (w  ,♦  ) at  t ■ 0. 

o o 

For  an  isotropic  radiator  (Lambertian  illumination)  at  t ■ 0 the  source 
function  becomes 
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B(t;u) 


P(t;|i,v1)  e T^W1 


(5) 


These  formulas  are  only  strictly  true  for  monochromatic  radiation, 
however,  with  suitable  averages  a wavelength  interval  may  be  simulated. 


The  boundry  conditions  are 

Ko;u,*)  - l(Tmax;-v,A)  - 0 

i.e.,  no  diffuse  radiation  incident  at  either  end  of  the  slab. 

Chandrasekkar  [1960]  showed  that  if  the  phase  function  could  be 
expressed  in  a Taylor-like-series  of  Legendre  polynomials  in  the 
scattering  angle,  then  by  the  addition  theorem  of  spherical  harmonics 
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the  phase  function  could  be  represented  in  the  Fourier  series  form 


N 

P(t;v,*;u V>  - I Pn(T;u,ti')  cos(n  - 1)  (*'  - *)  (6) 

n“l 

This  form  for  P suggests  that  I and  J could  be  similarily  expanded 
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I(t;v,+)  » Z In(T;M,u')  cos(n  - 1)  ($'  - $) 
n“l 


N 

» Z Jn(T ; u, u' ) cos(n  - 1)  (♦'  - ♦)  (7) 

n*l 


The  advantage  of  this  type  expansion  is  that  the  azimuthal  integral  in 
(2)  can  be  carried  out  analytically  with  considerable  reduction  in 
computation  time.  Furthermore,  since  the  forward  flux  at  t is  given  by 


Ff  (x ) 


I(T;u,$)udg  d$ 


(8) 


if  we  substitute  (7)  into  (8)  all  terms  with  n > 1 will  integrate  to 
zero.  Thus,  for  flux  calculations  we  need  only  consider  one  term  in  the 
Fourier  expansion  of  the  intensity  and  plase  functions. 


From  a distance,  a jet  engine  exhaust  has  characteristics  of  both  a 
mono-directional  source  and  a Lambert  source.  The  exhaust  gasses 
undoubtedly  emit  isotropically  but  at  a distance  only  those  photons 
which  come  from  a small  sterradian  will  be  detected.  For  this  reason  we 
have  made  calculations  with  both  sources  to  bracket  the  conditions 
existing  for  real  sources. 


Substituting  ♦'  ■ into  (7)  for  a monodirectional  source,  and  (7) 
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into  (1)  gives 


dl^x.y) 

dx 


.1 


where  J from  (2),  (4) 


I^x.y)  - w(x.)J*(x,y) 
and  (6)  is 


(9) 


Ax.y)  - 1/4  F e 


-r/u 


P (x;y,yQ)  + 


1/2/-, 


P^T, 


y »u 


') 


I 1 ( t , u ' )dy ' 

(10) 


Each  term  of  (9)  and  (10)  has  a cos  (n-1)  ($q-$)  multiplying  it  which 
can  be  neglected.  The  first  term  on  the  right  hand  side  of  (10) 
represents  the  contribution  to  the  source  function  due  to  primary 
scattering  of  the  direct  unidirectional  source  by  the  unit  volume  at 
level  x.  A successive  scattering  iteration  procedure  is  used  to  solve 
the  radiative  transfer  equation  (9)  using  this  first  term  as  the  initial 
value.  The  iteration  procedure  continues  until  values  of  I*(x,y)  on 
successive  iterations  converge  to  some  tolerance  for  all  values  of  |i. 
Numerically  (9)  is  solved  by  dividing  the  slab  into  a number  of  layers. 
Assuming  the  value  of  the  source  function  within  each  layer  to  be  the 
arithmetic  mean  of  the  values  at  the  edges  of  the  layer,  (9)  can  be 
written  1 

c 

ll(xl;+y)  - l1(''1+1;!»)e'ATk/M  - w(xk)  J1  (xk;+y)  [l-e~AVW] 

where  k ■ l for  +M  and  k * *— 1 for  -y.  I^(t^;u)  and  I*(x, ;-p)  represent 

t the  intensities  emergent  at  th^  front  of  layer  l-l  and  back  of  layer  t 

respectively.  All  values  of  I (x^ ;+y) , ^ forward  direction,  are  computed 
before  proceeding  to  computations  of  I (x,;-y).  ^This  means  a forward 
pass  from  the  source  outward  (increasing  t;  for  IA(x,y)  and. then  a 
backward  pass  (decreasing  P)  through  the  slab  calculating  I (x;-y). 

t The  integral  of  (10)  is  evaluated  by  the  trapezoidal  rule  of  integration 

using  increments  of  y corresponding  to  2°  increments  in  e (91  values  per 
layer) . 

C A similar  procedure  is  used  to  calculate  the  intensities  for  the 

isotropic  source  but  with  (5)  as  the  primary  source  function. 
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For  very  small  particles  (Rayleigh  scattering)  the  phase  function 
can  be  found  by  substituting  the  relationship  between  the  scattering 
angle  and  spherical  coordinates 

cos  e - uu*  ♦ ( 1 - y2)1/2  (1  - u,2>1/2  cos  ($'  - *) 
into  the  Rayleigh  scattering  expression  [Chandrasekkar , 1960] 

R(u,4;u',4')  ■ 3/4(1  ♦ cos2e) 
to  give 

R(u.u')  - 3/8  (3  - w2  - u'2  + 3iiV2) 

The  phase  function  for  large  particles  is  more  involved  and 
requires  much  more  computer  time.  Mie  [1908]  first  solved  for  the 
electromagnetic  field  scattered  by  a dielectric  sphere,  of  radius  and 
index  of  refraction  m = n^  - in2,  irradiated  by  a plane  wave.  The 
details  of  this  problem  in  electromagnetic  field  theory  are  nicely 
described  in  Chapter  13  of  Born  and  Wolf  [1970].  It  is  sufficient  to 
say  here  that  the  wave  intensity  scattered  into  a cone  of  unit 
solid  angle  vrtiose  axis  makes  an  angle  9 with  the  incident  beam  as 

X2 

I(x,m,e)  = ^ (sxz  (e)  + s22(6)) 

where 

OD 

Sl(e)  " z.  roTT  [Vn(cos  e)  + Vn(cos  9)1 

n*i 

00 

Sl(0)  “ £.  Ibnwn(cos  e)  + Vn(cos  0)1 

n*i 

and  *n  (cos  9)  and  rn  (cos  9)  are  expressed  in  terms  of  the  first  and 
second  derivatives  of  the  ordinary  Legendre  functions.  Compare  Chapter 
9 of  Van  De  Hulst  [1957]  and  Chu  and  Churchill  [1955],  Sekera  [1952] 
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and  Chu  and  Churchill  [1955)  shoved  that  the  scattering  functions  could 
be  simplified  by  repeated  use  of  recurrent  relationships  between  the 
derivatives  and  products  of  Legendre  functions.  This  would  yield  for 
these  functions  a series  whose  terms  are  ordinary  Legendre  polynomials 
weighted  by  coefficients  depending  upon  m and  x ■ 2 *a/A  only  [Dave, 
1970; 


A N 

I(x,m,e)  - S L^x.nOP^  (cos  0) 

Dave  J.V.  and  J.  Gazdag,  1970).  The  series  converges  for  all  values  of 
x and  m,  but  the  number  of  terms  required  depends  upon  the  size  of  x and 
m.  Generally  speaking  the  larger  x and  m,  the  greater  N.  The  Mie  phase 
function  is  then  calculated  by  generating  the  L^  for  specific  values  of 
x and  m.  The  legendre  functions  are  then  calculated  for  the  desired 
values  of  y and  y')  and  the  series  evaluated  to  find  M(y,y').  The  phase 
function  M(y,l)  is  plotted  in  figure  2 for  three  particle  sizes,  mie 
scattering  for  small  particles  is  seen  to  approach  a Rayleigh 
distribution  and  for  large  particles  is  primarily  a forward  scattering 
process.  Figure  3 illustrates  the  variation  of  the  scattering 
coefficient  (i.e.  the  ratio  of  t;otal  power  scattered  to  the  particle 
geometric  cross  section)  versus  x.  Similar  curves  are  found  in 
Deirmendjian  [1969). 

The  computer  code  will  actually  calculated  L^'s  for  a distribution 
of  particle  sizes.  The  coefficients  are  weighted  by  the  number  of 
particles  having  the  given  size  x.  For  the  calculations  presented  in 
this  report,  a gaussian  distribution  of  sizes  was  used  with  the  extreme 
radii  ♦_  5Z  of  the  center  radius. 

B.  Scattering  Particles  and  Absorbing  Gas 


t . 


The  first  task  was  the  calculation  of  the  radiation  flux  emerging 

from  a one  kilometer  thick  homogeneous  slab  of  pure  scattering 

particles.  The  "effective"  thickness  of  the  slab  depended  upon  the 

2 

particle  content  of  a one  cm  column.  The  transmission  through  the  slab 
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versus  column  content,  or  equivalently  density,  for  the  isotropic  and 
mon-directional  sources  are  plotted  in  Figure  4.  The  index  of 
refraction  for  these  calculations  was  m ■ 2.7  + iO,  which  corresponds  to 
the  real  part  of  the  carbon  refractive  index.  Actually  there  is  not 
much  difference  in  the  particle  scattering  effectiveness  for  m > 2.5, 
since  the  average  scattering  is  almost  independent  of  refractive  index 
above  about  m “ 2.5.  The  radiation  flux  at  the  source  was  equal  for 
both  cases.  Transmission  values  for  the  isotropic  source  were 
consistantly  less  than  the  monodirectional  source  because  the  average 
path  length  was  greater.  Since  a purely  scattering  medium  cannot 
destroy  photons,  Figure  4 implies  that  the  transmission  decreases  with 
increasing  optical  depth  because  increasingly  more  radiation  is 
backscattered . The  backscattered  radiation  is  the  difference  between 
100%  and  the  transmission,  so  that  the  sum  of  the  transmitted  and 
backscattered  intensities  is  unity.  This  is  just  the  statement  of 
conservation  of  energy. 

Based  on  these  characteristics  of  radiation  transport  in  a pure 

scatterer,  our  aim  was  to  evaluate  the  possible  synergistic  absorption 

in  a mixture  of  scattering  particles  and  absorbing  gas.  The  degree  of 

synergism  is  defined  as  the  amount  that  the  ratio  R =■  x /x  x is  less 

Pg  P g 

than  unity.  The  transmission  through  scattering  particles  alone  is  x , 

through  absorbing  gas  x , and  through  a mixture  of  scatterers  and 

absorbing  gas  x . (The  greek  symbol  x has  been  used  here  to  represent 

transmission  whereas  x was  used  in  the  computer  code  description  to 

represent  optical  depth.)  Scattering  particles  were  assumed  to  have  m * 

2.7  and  NjO  was  used  for  the  absorbing  gas.  Transmission  values  for 

three  particle  sizes,  x = .1,  1 and  10.  At  4.5um  this  corresponded  to  a 

particle  radium  of  .72,  7.2  and  72  microns.  Results  for  x ■ .1  are 

shown  in  Figure  5.  These  curves  illustrate  the  transmission  through  a 

one  km  thick  slab  for  three  values  of  scattering  particle  densities, 

corresponding  to  transmissions  of  20,  50  and  80%,  and  as  a function  of 

N2O  absorbing  gas  density.  The  solid  curves  represent  transmission 

through  1 km  of  scattering  x and  1 km  of  absorption  x , with  the  total 

p a 

transmission  x ■ x x . The  dashed  curves  in  Figure  5 are  the 
P a 

transmissions  through  a mixture  of  pure  scattering  particles  and 
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Fig.  4.  Transmission  through  a 1 km  thick  slab  of  pure  scattering  part 
both  uni-directional  and  isotropic  radiation  sources. 


absorbing  ^0  gas.  The  difference  between  these  two  sets  of  curves 
shows  that  some  synergism  does  occur.  On  an  absolute  scale  this  effect 
is  always  less  than  10%,  which  is  not  quite  enough  to  be  an  attractive 
infrared  obscuration  technique.  The  dot  dashed  curve  in  Figure  5 and 
the  synergism  ratio  R plotted  in  Figure  6 show  that  the  percent 
reduction  in  transmission  due  to  synergistic  effects  is  greatest  for 
small  values  of  transmission  (large  absorption).  However,  these 

8 -3 

relative  values  are  misleading.  If  we  take  a mixture  of  7.95  x 10  cm 

14  -3 

scattering  particles  and  8 x 10  cm  molecules  of  ^0  the  curves  show  a 

transmission  of  16%  through  a 1 km  thick  slab.  For  titanium  dioxide 

scattering  particles,  this  total  number  density  would  correspond  to  a 

“6  3 “6  3 

total  mass  density  of  5.06  x 10  gm/ cm  , made  up  of  5 x 10  gm/cm  for 

_8  3 13 

titanium  dioxide  and  6 x 10  gm/cm  for  N_0.  A density  of  1.8  x 10 

3 _ 2 

molecules/cm  of  N_0  would  also  give  a 16%  transmission  but  would  only 
-7^3 

weigh  1.3  x 10  gm/cm  . If  kerosene  had  the  same  scattering  properties 
it  would  only  reduce  the  scattering  particle  weight  by  a factor  of  4. 

We  conclude  from  these  calculations  that  pure  spherical  scattering 
particles  added  to  an  absorbing  gas  will  decrease  the  infrared 
transmission  through  the  gas,  but  unfortunately  the  effect  is  small  and 
unless  there  is  a very  light  weight  scattering  particle  the  combined 
weight  of  the  scattering/absorbing  mixture  will  exceed  that  of  just  an 
absorbing  gas.  Any  hope  to  use  this  technique  in  plume  IR  obscuration 
will  have  to  lie  with  scattering  particles  that  absorb  as  well  as 
scatter.  We  analyze  one  such  particle,  carbon,  in  the  next  section. 

C.  Carbon  Absorption 

The  complex  index  of  refraction  for  carbon  leads  to  both  absorption 
and  scattering  for  carbon  particles.  Since  the  value  of  the  refractive 
index  is  dependent  upon  the  physical  state  of  the  carbon  and  also  the 
wavelength  [Whitson,  1975]  we  have  used  m ■ 2.7  - il  which  is 
representative  of  soot  between  2.7  and  4.5  um.  For  a size  parameter  of 

-14 

x ■ .1,  the  volume  scattering  coefficient  is  only  B^  ■ 2.5  x 10  while 

the  volume  absorption  coefficient  is  over  1000  times  higher  at  8 * 1 x 

_1 1 t * 

10  . This  simply  means  that  carbon  particles  are  more  effective  as 

absorbers  than  as  scatterers.  A series  of  computer  runs  were  made  with 
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Fig.  6.  Synergistic  ratio  for  pure  scattering  particles 
and  N20  mixture  and  carbon-N20  mixture. 
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carbon  particles  of  sizes  x ■ .1,  1,  and  10  in  N^O  to  determine  possible 
synergistic  effects  due  to  scattering.  The  synergistic  ratio,  defined 
earlier,  for  some  of  these  calculations  are  plotted  in  Figure  6,  and 
show  that  the  carbon  scattering  does  not  increase  the  net  radiation  loss 
by  tfce  NjO  gas.  The  absorption  within  the  carbon  particle  is  too  great. 
This  absorption  is  illustrated  in  Figure  7 where  transmission  versus 
mass  density  for  three  different  particles  sizes  is  plotted.  Remember 
this  represents  transmission  through  a one-kilometer  slab  having  the 
uniform  mass  density  indicated.  It  is  most  interesting  to  see  that  x ■ 

1 is  a preferred  particle  size.  This  is  due  to  the  fact  that  the 
scattering  coefficient  attains  a maximum  value  near  x “ 1. 

The  effectiveness  of  carbon  particles  in  infrared  obscuration  is 
indicated  in  Figure  8.  IR  transmission  at  4.5Um  through  both  carbon  and 
NjO  is  plotted  versus  weight.  It  is  easily  seen  that  much  less  carbon 
is  required  for  a given  transmission  value  even  though  the  particles  are 
much  bigger  than  molecules. 

III.  EXPERIMENTAL  EFFORT 


A.  Experimental  Procedure 

The  apparatus  used  for  these  measurements  was  identical  to  that 
used  in  the  previous  work  [Harris  et  al.,  1975]  on  gaseous  absorbers, 
with  one  exception:  the  gas  absorption  cell  was  replaced  by  a 
scattering/absorbing  chamber.  This  chamber  was  a small  plexiglass  cube 
( 2"  inside  dimension)  fitted  with  saplire  windows  on  two  opposing 
faces.  A small  electric  motor  was  mounted  on  the  outside  of  the  top 
face . The  motor  shaft  extended  downward  through  the  box  and  drove  a 
multi-bladed  fan  fitted  to  the  end  of  the  motor  shaft.  Additionally,  a 
magnetic  stirring  bar  was  placed  on  the  floor  of  the  chamber  and  driven 
by  a magnetic  stirring  motor  from  the  bottom.  The  combination  of  fan 
and  stirring  bar  provided  a reliable  means  of  maintaining  reasonably 
uniform  particle  suspensions  within  the  chamber.  Neither  the  fan  nor 
stirring  bar  interfered  with  the  direct  optical  line  of  sight  through 


Transmission  through  1 km  slab  of  carbon  particles 


MASS(  1(5  GRAMS  /CM 


I 


the  chamber.  Early  attempts  to  suspend  the  particles  with  an  ultrasonic 
particle  disrupter  unit  did  not  meet  with  success.  The  worst  defect  of 
this  method  was  the  generation  of  standing  wave  patterns  in  the  chamber 
with  the  result  of  highly  non-uniform  particle  distributions.  In  one 
columnar  chamber  design,  for  example,  the  particles  noticibly  collected 
at  the  wave  nodes.)  No  quantitative  measure  of  particle  number  density 
was  made.  This  parameter  is  eliminated  in  the  method  of  analysis  used 
(described  later)  and  is  not  required.  Crude  estimates  of  maximum 
density  obtainable  in  the  chamber  were  made  by  noting  visually  the 
largest  amount  of  material  that  could  be  maintained  off  the  floor  of  the 
chamber.  This  maximum  was  approximately  0.00 l g/cm  . Two  copper  tubes 
through  the  top  of  the  chamber  were  used  to  flush  in  ^0  and  N2  gases . 
The  radiation  source,  chopper,  atmospheric  simulation,  monochrometer, 
detector  and  display  apparatus  were  the  same  as  used  in  the  previous 
work,  [Harris  et  al.,  1975]. 


Particles  used  in  the  measurements  included  carbon  in  the  size 
category  0.1  to  0.66  um  diameter  (0.32  urn  average,  x * .22  at  4.5pm)  and 
3 to  8 pm  diameter  (6  pm  average),  and  2 pm  diameter  TiC^* 

Data  were  taken  for  three  conditions  in  the  scattering  chamber: 


1)  chamber  empty  (flushed  with  N2) 

2)  NjO  alone 

3)  particles  in  N20. 

Typical  spectral  results  in  the  4.6  um  region  are  shown  in  Figures  9 and 
10.  The  data  were  analyzed  on  a whole-band  basis  by  measuring  areas 
under  the  various  curves.  Measurements  for  conditions  1 and  2 provide 
the  transmission  data  for  N^O  alone.  Measurements  for  conditions  2 and 
3 provide,  simultaneously,  the  transmission  data  for  scatterers  alone 
and  scatterers  and  absorption  gas  in  combination.  This  was  accomplished 
by  determining  the  transmissions  separately  in  the  4.2  um  spike  and  the 
4.3  to  4.7  um  wing.  Since  there  is  no  absorption  due  to  NjO  in  the  4.2 
um  region,  measurements  on  this  spike  provide  a measure  of  attenuation 
due  to  scatterers  alone.  This  measure  of  particle  attenuation  is 
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Fig.  10.  Upper  spectrum,  solid  curve,  is  4.1  to  4.8  micron  region  of 
flame  as  seen  through  C0_  atmospheric  path  simulator  (sea 
level,  1200ft.)  with  600  torr  partial  pressure  of  N_0  in  scat- 
tering cell.  Dashed  curves  are  for  two  different  concentrations 
^ of  carbon  in  the  N„0  mixture. 

Lower  spectrum  is  tor  the  same  region  and  same  path  length 
simulation  with  nitrogen  in  scattering  cell  (solid  curve)  and 
nitrogen  plus  two  different  concentrations  of  suspended  carbon 
(dashed  curves). 
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assumed  to  prevail  throughout  the  entire  band  spectral  region. 
(Preliminary  measurements  of  particle  attenuation  had  been  made  by 
suspending  the  particles  in  Nj . It  was  discovered,  however,  that  for 
fixed  stirrer  power  levels,  that  Nj  and  NjO  would  not  suspend  the  same 
amount  of  material.)  Measurements  on  the  4.3  to  4.7  Km  wing  provide  the 
measure  of  attenuation  due  to  particles  and  NjO  in  combination. 

Let  Tg  be  the  measured  transmission  of  the  NjO  alone  (mean  value 
over  whole  band),  tp  be  the  measured  transmission  due  to  particles  alone 
(mean  value  in  4.2  pm  spike),  and  tpg  be  the  measured  transmission  due 
to  particles  and  NjO  in  combination  (mean  value  in  4.3  to  4.7  pm  wing). 

B.  Experimental  Results 

Results  for  the  synergism  ratio  R are  presented  in  Figures  11  and 
12  for  the  0.3  pm  carbon  sample  and  the  2 pm  Ti02  sample,  respectively. 
The  plotting  abscissa  for  these  figures  is  the  particle-alone 
transmittance  Tp.  To  within  the  experimental  uncertainty,  the  ratio  R 
is  unity  for  all  optical  depths.  The  scatter  of  the  data  is  quite  large 
and  could  conceivably  hide  a synergism  of  less  than  R “ 0.9  for  carbon 
(none  for  Ti02).  However,  this  amount  of  synergism  would  not  be 
significant  in  the  overall  plume  obscuration  scheme. 

Data  obtained  in  the  2.9  um  region  and  data  for  other  particle  size 
groups  were  too  scattered  and  incomplete  for  presentation.  Given  the 
negative  result  of  Figures  11  and  12,  it  was  felt  that  no  purpose  would 
be  served  be  spending  time  on  improving  that  data. 
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Fig.  11.  Synergism  ratio  for  0.3pm  carbon  sample 
in  the  4.3pm  spectral  region. 


Fig.  12.  Synergism  ratio  for  2.0pm  TiC^  sample 
in  the  4.3pm  spectral  region. 


IV.  CONCLUSIONS 
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A study  was  conducted  to  evaluate  the  role  of  scattering  particles 
in  the  absorption  of  infrared  radiation  through  an  aerosol-absorbing  gas 
mixture.  The  study  incorporated  both  theoretical  and  experimental  methods. 
Theoretical  calculations  of  IR  radiation  transport  through  scattering/ 
absorbing  mediums  were  made  with  a computer  code  obtained  from  Lawrence 
Radiation  Laboratories  and  modified  by  us.  Experimental  measurements  were 
made  of  IR  flame  spectra  transmitted  through  cells  filled  with  scattering 
particles  and/or  absorbing  gases.  The  results  of  these  two  efforts  are 
consistent  and  may  be  summarized  as  follows: 
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1)  . Both  theoretical  and  experimental  results  show  that 

synergistic  effects  in  media  containing  pure  scatterers 
and  absorbing  gas  is  small,  not  greater  than  about  10% 
Titanium  dioxide  was  used  for  the  scattering  medium  and 
^0  the  absorbing  gas  because  of  its  absorption  effective- 
ness in  the  red  wings  of  the  2.7  and  4.3ym  bands. 

2) .  Theoretical  and  experimental  results  also  showed  that 

synergistic  effects  in  a medium  containing  carbon  and  N20 
was  again  small,  less  then  10%.  Carbon  particles  both 
scatter  and  absorb  internally  infrared  radiation. 

3) .  In  the  light  of  conclusions  1 and  2,  we  regretfully 

conclude  that  augmented  absorption  produced  by  scatter- 
ing particles  in  an  absorbing  medium  is  not  a viable 
technique  for  infrared  observation  in  exhaust  plumes. 

4) .  Carbon  particles  produce  strong  infrared  absorption. 

The  internal  absorption  coefficient  is  about  10 J higher 
than  the  scattering  coefficient.  Because  the  absorption 
is  particle-like,  there  is  less  structure  in  the  absorp- 
tion-versus-frequency  response  than  in  a gaseous  absorber 
like  N-0.  Carbon  absorption  at  2.7ym  is  only  a factor 
of  2 less  than  that  at  4.3ym.  Furthermore,  it  was  shown  at 
4.5ym  that  a greater  weight  of  ^0  was  required  to  achieve 
a given  IR  absorption  than  with  carbon.  Except  for  the 
black  color  in  the  visable,  black  carbon  is  an  attractive 
substance  for  hot  plume  obscuration. 
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